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The role of shear in alleviating catastrophic quenching by shedding small-scale magnetic helicity through fluxes along 
contours of constant shear is discussed. The level of quenching of the dynamo effect depends on the quenched value of 
the turbulent magnetic diffusivity. Earlier estimates that might have suffered from the force-free degeneracy of Beltrami 
fields are now confirmed for shear flows where this degeneracy is lifted. For a dynamo that is saturated near equipartition 
field strength those estimates result in a 5-fold decrease of the magnetic diffusivity as the magnetic Reynolds number 
based on the wavenumber of the energy-carrying eddies is increased from 2 to 600. Finally, the role of shear in driving 
turbulence and large-scale fields by the magneto-rotational instability is emphasized. New simulations are presented and 
the 37r/4 phase shift between poloidal and toroidal fields is confirmed. It is suggested that this phase shift might be a 
useful diagnostic tool in identifying mean-field dynamo action in simulations and to distinguish this from other scenarios 
invoking magnetic buoyancy as a means to explain migration away from the midplane. 
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1 Introduction 

Shear clearly plays an important role in amplifying toroidal 
fields from poloidal, but that is not all. Shear also plays a 
role in "unquenching" any dynamo effect that may play a 
role in producing poloidal field from toroidal, thus closing 
the dynamo loop. A prime example of such a dynamo effect 
is the a effect, but other possible known effects may include 
the shear-current effect and the incoherent alpha-shear ef- 
fect. The "unquenching" of dynamo effects, as well as the 
dynamo effects themselves, require more detailed consid- 
erations of the results available so far. This is the principal 
goal of this paper. 

The term "unquenching" in connection with the a ef- 
fect may appear somewhat unusual, but this choice of words 
must be seen in contrast to the possibility of catastrophic a 
quenching. Here, "catastrophic" indicates that the quench- 
ing by the mean magnetic field, B, becomes more extreme 
as the magnetic Reynolds number, R m , increases. Tradition- 
ally, such quenching is represented by the simplistic for- 
mula (Vainshtein & Cattaneo 1992) 



a(B) 



1 + R m B /B 2 q 



(1) 



where B cq — (nopu 2 ) 1 / 2 is the equipartition field strength 
with respect to the kinetic energy density, u is the small- 
scale turbulent velocity, p is the density, po is the vacuum 
permeability, R m = u rms /rjk{ is the magnetic Reynolds 
number, rj is the magnetic diffusivity, k{ is the wavenum- 
ber of the energy-carrying scale, and angular brackets de- 
note volume averaging. However, Eq. (1) is really only valid 



under special circumstances that are quite uninteresting for 
dynamo action: infinite wavelength of the magnetic field, 
complete stationarity, and no possibility of magnetic helic- 
ity fluxes. The latter possibility is now believed to be the 
most important one for astrophysical dynamos, as was first 
suggested by Blackman & Field (2000). Any one of these 
three caveats alleviates the severity of catastrophic quench- 
ing, because they all lead to "extra terms" that enter in the 
numerator of Eq. (1) with an R m factor in front, i.e. 
&o + Rm x "extra effects" 



a(B) = 



(2) 



1 + R m B /Bg q 
where we have assumed that the kinematic a value, ao, re- 
mains independent of time. In its essence, this equation with 
extra effects included goes back to the early work of Klee- 
orin & Ruzmaikin (1982), and later Kleeorin et al. (1995, 
2000), and has been discussed in connection with allevi- 
ating catastrophic a quenching by Blackman & Branden- 
burg (2002) and Brandenburg & Subramanian (2005a). A 
more complete quenching formula with extra effects in- 
cluded takes the form 



a + R B 



HqJ-B 

*h b'L 



da /dt 
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(3) 

1 + R m B" /B 2 q 
Let us now discuss separately all three terms in the paren- 
thesis of the numerator of Eq. (3). 

(i) The effect of the J ■ B term is clearly seen when 
considering the saturation of homogeneous dynamos in a 
periodic domain. In that case this formula gives 

a(B) -> rjtpoJ ■ B/B 2 (fori? m ^oo), (4) 
so there is nothing catastrophic about this formula, unless 
rjt itself is catastrophically quenched. Of course, if mean 
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fields are defined as full volume averages, B becomes com- 
pletely uniform, so /j,qJ = V x B = 0. This is the case in 
numerical experiments by Cattaneo & Hughes (1996). The 
other case with a finite J ■ B term was seen in the simula- 
tions of Brandenburg (2001), where not only a large-scale 
magnetic field was found to saturate at super-equipartition 
values, but also the a and r)t effects were found to be only 
mildly quenched. 

(ii) The d/dt term in Eq. (3) is important to explain the 
absence of an otherwise premature onset of quenching at re- 
sistively low field strengths where B /B^ q « i?~ 1 . This is 
also confirmed by controlled numerical experiments where 
the initial field was a weak Beltrami field (Brandenburg et 
al. 2003). 

(iii) Finally, the effect of magnetic helicity fluxes was 
first seen in simulations with imposed fields by Branden- 
burg & Sandin (2004), and then later in dynamo simulations 
(Brandenburg 2005a), which brings us to the main topic of 
this paper. As was already seen in earlier simulations with- 
out shear, just allowing for open boundary conditions alone 
does not help to produce a finite magnetic helicity flux in 
Eq. (3) and hence does not alleviate catastrophic a quench- 
ing (Brandenburg & Dobler 2001). However, Vishniac & 
Cho (2001) showed that in the presence of differential ro- 
tation a magnetic helicity flux can be generated and that it 
flows along the rotation axis. More detailed work of Sub- 
ramanian & Brandenburg (2004, 2006) and Brandenburg & 
Subramanian (2005b) resulted in a simple formula for this 
particular contribution to the flux: 

Tc = C vc (SB)xB, (5) 

where = \{Ui.j + Uj^) is the rate of strain matrix of 
the mean flow, and Cyc is a dimensionless number of order 
unity. This flux is along contours of constant shear, as was 
demonstrated by Brandenburg et al. (2005); see also Fig. 1. 

The sign of the magnetic helicity flux is negative in 
the Northern hemisphere, so vectors of positive magnetic 
helicity flux point away from the surface in the Northern 
hemisphere, then through the equator and into the Southern 
hemisphere. 

2 Quantitative considerations 

Simulations of Brandenburg & Sandin (2004) indicate that 
the magnitude of the Vishniac-Cho flux might be on the or- 
der of 10 46 Mx 2 /cycle, if applied to the Sun. In units of 
To = u lms k[ Bq the nondimensional flux was estimated to 
be about 30/^Sh, where 

Sh = S/u rms k{ (6) 

is the nondimensional shear parameter and S is the shear 
rate. The corresponding Cyc parameter would then be 1 . 
This result appears in conflict with the expectation that 
Cyc should be of order unity (Subramanian & Branden- 
burg 2004; Brandenburg & Subramanian 2005b). It would 
therefore be important to return to this issue using a simpler 




Fig. 1 Vectors of Tq together with contours of U y which also 
coincide with the streamlines of the mean vorticity field W. Note 
the close agreement between Tc vectors and W contours. The 
orientation of the vectors indicates that negative current helicity 
leaves the system at the outer surface (a; = 0). Adapted from Bran- 
denburg et al. (2005). 

one-dimensional shear profile together with open boundary 
conditions. 

A more conclusive indication for the operation of helici- 
ty fluxes is the demonstration of a successful dynamo sim- 
ulation; see, e.g., Brandenburg (2005a), where it was also 
demonstrated that with closed perfect conductor bound- 
ary conditions the dynamo was not successful in produc- 
ing large-scale fields (at least not during the course of the 
present simulation). In Fig. 2 we plot the ratio of the en- 
ergy contained in the large-scale field to the total magnetic 
energy. We see clearly that large-scale dynamo action is 
only possible with open boundary conditions. By open we 
mean here the vertical field condition that is commonly used 
in simulations of magneto-convection (Hurlburt & Toomre 
1988). Such boundary conditions do permit a finite mag- 
netic helicity flux, but they do not allow Poynting flux to 
pass through the boundaries. It is at present unclear whether 
the absence of a Poynting flux is a serious short-coming or 
not. 

Another remarkable result found in the work of Bran- 
denburg (2005a) is that large-scale magnetic fields can be 
found even without kinetic helicity and hence without a ef- 
fect. An obvious possibility might be that this is caused by 
the shear-current effect of Rogachevskii & Kleeorin (2003, 
2004). However, in the case of a simpler geometry where 
the contours of constant shear are purely vertical, which is 
relevant for accretion discs, for example, no direct support 
for the existence of this effect has been found (Branden- 
burg 2005b; Brandenburg et al. 2008a). This is also consis- 
tent with earlier analytic results of Riidiger & Kitchatinov 
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Fig. 2 Evolution of the ratio of the energies contained in the 
large-scale field to the total magnetic energy for open and closed 
boundaries. Note that large-scale dynamo action is only possi- 
ble with open (vertical field) boundary conditions. Adapted from 
Brandenburg (2005a). 

(2006) and Radler & Stepanov (2006) using the second or- 
der correlation approximation (SOCA). In essence, a shear- 
current effect would be described by a mean-field equation 
with an anisotropic magnetic diffusion tensor, r]ij . We write 
the governing equation here for the magnetic vector poten- 
tial, 

Y)~A — — 

-j^- = -AjUj,i - fJLOTjijJj - (MtrjJi, (7) 

where J = — V 2 A, and one-dimensional averages have 
been employed, i.e. J = J(z, t) in the present case, and 
U = (0, Sx, 0). In the case of vertical contours of the mean 
shear, we have C/2,1 = S and all other components van- 
ish. Here, the cross-stream direction is i = x or 1, and the 
streamwise direction is i = y or 2. 

In order to have a closed dynamo loop, one would then 
need to have a finite 7721 component with the same sign as 
that of S. According to SOCA calculations and simulations 
this is however not the case. In Fig. 3 we show the resulting 
values of 7721 for such a linear shear flow, as obtained using 
the testfield method. 

Here we have used S < 0, so we are looking for neg- 
ative values of 7721 for an operational shear-current effect. 
In all cases we find 7721 > 0, except for large values of R m 
when P m — 20. However, the error bars are large. There is 
perhaps the possibility that the sign of 7721 may change un- 
der other circumstances, e.g. for the more complicated shear 
profile shown in Fig. 1. Yet another possibility is that in the 
presence of helicity the sign may change. Some evidence 
to this effect has been provided by Mitra et al. (2008). It 
turns out that in nonhelical shear flow turbulence, 7721 can 
become negative in the saturated state with helicity. How- 
ever, given that there is still an a effect, the relevance of the 
shear-current effect is less obvious in these simulations. 

Obviously, with helicity there is also an a effect, so the 
shear-current effect would not be the sole cause of large- 
scale dynamo action. This leaves us with the question what 
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Fig. 3 Dependence of 7721, normalized by the total (turbulent 
and microscopic) magnetic diffusivity on shear parameter (upper 
panel) and on magnetic Reynolds number for Re = 1.4 (middle 
panel) and for a fixed magnetic Prandtl number, P m = 20 (lower 
panel). Adapted from Brandenburg et al. (2008a). 

causes large-scale magnetic field generation in non-helical 
turbulence with shear? In addition to the simulations of 
Brandenburg (2005a), discussed above, there are also simu- 
lations of Yousef et al. (2008), where a large-scale magnetic 
field is generated. They find a growth rate that is propor- 
tional to the shear rate S. The authors argue that such a re- 
sult would not be consistent with the shear-current effect, 
because the growth rate would then be proportional to the 
product of S and 7721, where 7721 itself would be proportional 
to 5. 

The alternative proposal by Brandenburg et al. (2008a) 
is that an incoherent a-shear or aQ dynamo is at work. 
The occurrence of an incoherent a-shear dynamo (cf. Vish- 
niac & Brandenburg 1997; Proctor 2007; Kleeorin & Ro- 
gachevskii 2008) was quantified by estimating the rms val- 
ues of the fluctuations of all components of a>ij and r]ij. 
Figure 4 shows how these rms values vary with increasing 
values of R m . We recall that the average of all components 
of ctij is zero, so there is no regular a effect. The onset of 
incoherent dynamo action depends on the value of the dy- 
namo number 

^lT h = a rms S/ VT kl (8) 
where a rms is the rms value of the streamwise component 
a22 (but all components are found to have the same rms 
value), t]t = f]T + V is the sum of turbulent and micro- 
scopic magnetic diffusivity, and k\ — 2n/L is the lowest 
wavenumber in the z direction of the box. 
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Fig. 4 Dependences of the rms values of the temporal fluctua- 
tions a rms (normalized by r/tokt), r]l ma , 77™ s , and rf™ (normal- 
ized by ?jto), on R m for Re = 1.4 and Sh = —0.6. Adapted from 
Brandenburg et al. (2008a). 



The simulations of Brandenburg et al. (2008a) give fluc- 
tuations of a rms that correspond to D™g oh « 4. This is to 
be compared with the marginal value of w 2.3 obtained nu- 
merically using a single-mode approximation. There is in 
principle also the incoherent shear-current effect, based on 
the random occurrence of values of 7721 with suitable sign. 
However, the effect is found to be subdominant compared 
with the dynamo number of the incoherent a-shear dynamo. 
Yet another possibility for generating large-scale magnetic 
fields in the presence of shear is given by additive noise 
in the electromotive force due to small-scale dynamo ac- 
tion (Blackman 1998). However, in Yousef et al. (2008) no 
small-scale dynamo was excited, so this type of noise would 
not explain the field seen in their simulation. 

In all cases one would expect random reversals of the 
toroidal magnetic field on a turbulent diffusion time scale. 
However, simulations usually give longer time scales, which 
may be explained by a tendency to conserve magnetic helic- 
ity. Indeed, using magnetic helicity conservation, Branden- 
burg et al. (2008a) found that in a closed domain the reversal 
time increases with i? m to the 1/2 power. The mean reversal 
time might therefore well be much larger than a turbulent 
diffusion time. 



3 Shear in helical turbulent dynamos 

Let us now look at the effects of shear in simulations 
where kinetic helicity and thus an a effect are present. 
In these cases there is the possibility of oscillatory solu- 
tions with frequency w cyc and propagating dynamo waves 
of wavenumber k and wave speed c = uj cyc /k. The sense 
of propagation is determined by the sign of the product aS. 
When shear is sufficiently strong (\S/ak\ 3> 1 or, equiv- 
alently, |Sh| 3> gfci/fcf for fully helical turbulence with 



a 



3 w rmsX the magnitude of the oscillation frequency is 



given by (akS/2) 1 / 2 , but in the marginal state this must 
be balanced by the diffusion rate r/xfc 2 - As was stressed 
by Blackman & Brandenburg (2002), this provides there- 
fore a robust tool for determining empirically the quenched 
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Fig. 5 Rm dependence of rjt in the dynamo-saturated state as 
measured by the testfield method (solid line; Brandenburg et al. 
2008b), and estimated for oscillatory aO, dynamos in a shearing 
box (dashed line; Kapyla & Brandenburg, in preparation). 



value of t]t = rjt + r\. However, it has been rather hard to 
reach large enough values of the ratio rjt jf] (which is also a 
proxy of i? m ) to make conclusive statements about the In- 
dependence of r] t . Indeed, the values obtained so far (Fig. 5) 
confirm recent results of Brandenburg et al. (2008b), where 
a measure of rj t has been obtained using the testfield method 
applied to the nonlinear state for i? m ranging between 2 and 
600. Unfortunately, those results leave some ambiguity be- 
tween the a effect and the turbulent magnetic diffusivity, 
because of terms of the form (J • B)B, whose coefficients 
can be associated both with a and with rjt. 

4 MRI-driven turbulence 

Interesting systems where the effects of shear play an ab- 
solutely vital role are seen are local simulations of accre- 
tion disc turbulence that is driven by the magneto-rotational 
instability (Hawley et al. 1995). When there is also strati- 
fication about the midplane, the Coriolis force produces an 
a effect. Simulations of Brandenburg et al. (1995) give rise 
to large-scale dynamo action that turns out to be oscillatory 
with dynamo waves propagating away from the midplane. 
In Fig. 6 we show new simulations that are similar to those 
of Brandenburg et al. (1995), but with somewhat larger res- 
olution and without using hyperviscosity. Only shock vis- 
cosity and shock resistivity are being used. Another dif- 
ference is that these simulations also have a potential field 
boundary condition instead of the vertical field condition 
that was used in Brandenburg et al. (1995). Given the re- 
markable similarity with the earlier simulations, we must 
conclude that the different boundary conditions do not seem 
to have a major effect. The new simulations are similar in 
that they too show oscillations with a typical period of about 
30-50 orbits (the orbital time is defined as T rot = 2ir/il). 
However, the new simulations also show considerably more 
fluctuations with parity variations and, more importantly, 
a noticeable decoupling of behavior in the Northern and 
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Fig. 6 (online colour at: www.an-journal.org) Space-time (z- 
t) diagram of B v (z,t for a run where turbulence is driven by 
the magneto-rotational instability. The magnetic Reynolds num- 
ber, based on the actual rms velocity and the pseudo wavenumber 
k f = 2-k/H, is about 120. 
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Fig. 7 (online colour at: www.an-journal.org) Mean poloidal and 
toroidal fields in a narrow slice near the upper disc plane as a func- 
tion of t and z, together with the averages over the same slice 
showing the phase relation between the two. The solid line indi- 
cates B y /B cq while the dotted line indicates B x /B e q, scaled up 
by a factor of 20 to make it better visible. 

Southern disc planes. The parity varies therefore between 
more nearly symmetric (even) and more nearly antisymmet- 
ric (odd) parity. 

5 Detecting mean-field dynamo action 

A possible means of identifying ail-type dynamo action as 
the main course of oscillations seen in simulations we pro- 
pose to determine the phase relation between poloidal and 
toroidal fields. This used to be a standard tool in solar dy- 
namo theory to infer the sense of radial differential rotation, 
but may also become an important tool in disc and other os- 
cillatory dynamos. Mean-field theory predicts a phase shift 
by jTT, which was confirmed by Brandenburg & Sokoloff 
(2002). Another alternative explanation for the migration 




-0.04 -0.02 0.00 0.02 0.04 

Fig. 8 Phase plot of the averages of poloidal and toroidal fields 
over the same narrow slices as in Fig. 7. The first times are plot- 
ted as dashed lines, but the later times are solid with increas- 
ing thickness toward the end showing that a point on the curve 
moves forward in a clockwise direction. Overplotted are two el- 
lipses showing B y o oc cos(ujt + (j>) versus B x o oc cos cut with 
4> = 0.65tt, B x0 = 0.035 B cq , and B x0 = 1.8 B cq (solid line) 
and (f> = 0.757T, B x0 = 0.03 B cq , and B x0 = 1.74 B oq (dotted 
line). 

away from the midplane would be magnetic buoyancy. This 
was discussed by Vishniac & Brandenburg (1997), who 
noted that the migration speed is only about 3% of the tur- 
bulent rms velocity. The idea of magnetic buoyancy playing 
a leading role was expressed again in connection with re- 
cent simulations of Blaes et al. (2008), where the outer disc 
surface is marked by a surface where the optical depth in a 
radiative transfer calculation was of order unity. However, 
no detailed proposal for the phase relation from the buoy- 
ancy effect has yet been made. 

Brandenburg et al. (1995) found that the details of the 
large-scale magnetic field generation can be described by an 
a effect that is negative in the Northern disc plane. Newer 
determinations of ay and T)ij using the testfield method 
confirm this result also for the present simulations; see 
Fig. 9. Again, we have here access to all 4+4 components 
of the ctij and rjij tensors. As in Brandenburg & Sokoloff 
(2002), who used a correlation method instead of the test- 
field method, a yy is negative in the upper midplane and has 
an extremum atzw ±H . However, an (= a xx ) is positive 
in the upper disc plane. 

In all cases we find that rjn « 7722 and always posi- 
tive, in contrast to earlier work using the correlation method. 
This function has a minimum in the midplane and grows 
away from the midplane by a factor of about 4. The off- 
diagonal components of and rjij are shown in Fig. 10. 

It turns out that, like a xx and a vy , also a xy and a yx are 
approximately antisymmetric about the midplane and posi- 
tive in the upper disc disc plane. The off-diagonal compo- 
nents of ay are normally interpreted in terms of turbulent 
pumping with an effective vertical velocity 7 = 5 (oi yx — 
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Fig. 9 Vertical profiles of a xx , ot yy (first and second panel), both 
normalized by the turbulent rms velocity, as well as r\ xx (solid line) 
and rjyy (dashed line) for an MRI disc simulations. 




Fig. 10 Vertical profiles of a xy , ot yx (first and second panel), 
both normalized by the turbulent rms velocity, as well as ij xy (solid 
line) and rj yx (dashed line) for an MRI disc simulations. 



a xy ). However, since a yx and a xy have the same sign, there 
is some cancelation. Nevertheless, since a yx is larger than 
a xy , there is net transport away from the midplane. This 
concerns predominantly the toroidal field component, B y , 
while the poloidal component, B x , is transported predom- 
inantly toward the midplane. Such differential pumping of 
poloidal and toroidal fields was first discussed by Kitchati- 
nov (1991) and later confirmed in simulations by Ossendri- 
jver et al. (2002). 

The off-diagonal components of rjij are symmetric 
about the midplane and both are mainly positive. The r] yx 
component is important for the shear-current effect, but it is 
found to have the wrong sign for this effect to operate in our 
simulations; see also Sect. 2. 

6 Dynamos in shearing convection 

Large-scale dynamo action of significant amplitude has 
been seen in global simulations of the geodynamo (Glatz- 
maier & Roberts 1995). However, simulations in Cartesian 
boxes have not yet shown such behavior. A recent exam- 
ple is that by Tobias et al. (2008) where strong shear was 
present. However, the absence of large-scale fields of signif- 
icant magnitude may be related to the orientation of shear. 
In their case those contours are horizontal, so the Vishniac- 
Cho flux, that goes along those contours, would not be able 
to escape the horizontally periodic domain. If it is true that 
large-scale dynamo action only works efficiently if excess 
small-scale magnetic helicity is expelled through the bound- 



aries, then it would be plausible that shear may not help the 
dynamo in all cases. More recent convection simulations by 
Kapyla et al. (2008) do show large-scale magnetic field in 
the presence of shear. They have contours of constant shear 
that do cross the boundaries. In fact, their shear is just like 
that in the accretion disc simulations discussed earlier. 

7 Discussion 

In this paper we have illustrated the important dual role 
played by shear: helping to produce strong toroidal field and 
helping to unquench the a effect. We have tried to make the 
case that successful large-scale dynamos must have an op- 
portunity to shed small-scale magnetic helicity through the 
boundaries. This is what the Sun does (e.g., Demoulin et al. 
2002), and this is also what many simulations can do, al- 
though the degree of realism to which simulations can do 
this varies. So far, there is no simulation that goes all the 
way to including at least a simplified way of modeling coro- 
nal mass ejections as the final stage in shedding small-scale 
magnetic helicity. This should obviously be tried in the near 
future using simulations in spherical shell segments, possi- 
bly with an additional outer layer where wind acceleration 
can occur (Brandenburg et al. 2007). 

In this connection it might be helpful to explain the spe- 
cial meaning of "small scale" in connection with large-scale 
dynamos. In solar physics, one is used to associating active 
regions with large scales, while small scales would often 
refer to the resolution limit of 100 km or less. Obviously, 
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in connection with the solar cycle, relevant time and length 
scales would be on the order of years and several hundred 
megameters, respectively - well encompassing the duration 
and size of active regions. In the present work we have used 
averages that do not a priori separate between large and 
small scales. By assuming an average over one or two coor- 
dinate directions (e.g. an azimuthal average for global sim- 
ulations of the Sun, or horizontal averages in local box sim- 
ulations), there could still be residual fluctuations on short 
time and small length scales, but they should go away with 
increasing system size in the azimuthal or horizontal direc- 
tions like one over the square root of the number of eddies 
that are being averaged over. Exactly this fact is behind the 
operation of the incoherent a-shear dynamo discussed at the 
end of Sect. 2. Clearly a finite horizontal extent is not only 
a numerical restriction of local boxes. Indeed, a finite az- 
imuthal extent is bare reality even for the Sun. This point 
has been made by Hoyng (1993) in order to explain limits 
to the phase and amplitude stability of the solar cycle, but 
this very mechanism could even suffice to drive a dynamo 
in shearing systems where there is no a effect present. It 
should be emphasized that the incoherent a effect as such 
is independent of shear, but it is able to produce coherent 
large-scale fields only in the presence of shear. Obviously, 
more work is needed to sharpen the analytical and numer- 
ical tools to understand such systems better. Finally, let us 
emphasize that incoherent a-effect dynamos too have a non- 
linear stage that is controlled by magnetic helicity at some 
level. In Brandenburg et al. (2008a) we only used a very 
simplistic one-mode truncation to make this point, but the 
whole dynamical quenching formalism applies otherwise 
just as well. 
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